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LONG-TERM  GOAL 

The  overall  goal  of  this  work  is  to  take  oceanic  radiative  transfer  theory  into  a  new  domain:  optically 
shallow  waters  with  spatially  variable  bottom  optical  properties,  including  sloping  bottoms  and  non- 
Lambertian  bottom  bidirectional  reflectance  distribution  functions. 

OBJECTIVES 

Currently  available  models  for  analysis  of  hyperspectral  remote  sensing  imagery  in  shallow  waters 
almost  always  assume  that  the  bottom  is  a  level,  homogeneous,  Lambertian  reflecting  surface — a 
spatially  homogeneous  surface  at  a  constant  depth  whose  reflected  radiance  appears  the  same  from  all 
viewing  directions.  The  first  objective  of  this  year’s  work  was  to  quantify  the  extent  to  which  non- 
Lambertian  bottoms  can  affect  upwelling  radiances  as  would  be  detected  by  above- water  sensors. 
Many  studies  of  marine  light  fields  also  assume  the  water  and  bottom  optical  properties  to  be 
horizontally  homogeneous,  which  is  not  the  case  in  shallow  waters  with  patchy  or  sloping  bottoms. 

The  second  objective  this  year  was  to  quantify  the  extent  to  which  patchy  or  sloping  bottoms  influence 
the  water-leaving  radiance.  In  particular,  I  wished  to  see  when  a  1-D  (horizontally  homogeneous, 
constant  depth)  radiative  transfer  model  is  sufficient,  and  when  a  computationally  expensive  3-D 
Monte  Carlo  model  must  be  used. 

APPROACH 

The  Hydrolight  radiative  transfer  model  (http://www.sequoiasci.com;  Mobley  and  Sundman,  2001a,b) 
allows  for  non-Lambertian  bottom  boundaries.  However,  because  of  the  lack  of  measured 
bidirectional  reflectance  distribution  functions  (BRDFs;  Mobley,  Zhang,  and  Voss,  2001;  Zhang  et  ah, 
2001)  for  actual  ocean  bottom  materials.  Hydrolight’s  mathematical  capability  had  not  been  previously 
exploited.  This  year,  non-Lambertian  BRDFs  were  used  in  Hydrolight  to  simulate  the  effects  of 
horizontally  uniform  but  non-Lambertian  bottoms  on  in- water  and  water- leaving  radiances. 

The  Hydrolight  model  requires  the  water  and  the  surface  and  bottom  boundaries  to  be  horizontally 
homogeneous;  Hydrolight  is  therefore  not  applicable  to  shallow  waters  with  spatially  variable  or 
sloping  bottoms.  To  simulate  the  inherently  3-D  light  fields  due  to  the  effects  of  spatially  variable  or 
sloping  bottoms,  Monte  Carlo  models  must  be  used. 
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WORK  COMPLETED 


This  last  year  of  CoBOP  work  was  used  to  conclude  the  model  development  and  studies  begun  in 
earlier  years,  and  to  prepare  and  submit  papers  on  that  work.  In  particular,  non-Lambertian  BRDFs 
were  used  to  simulate  radiances  for  situations  of  interest  to  the  airborne  Portable  Hyperspectral 
Imaging  Low-Light  Spectrometer  (PHILLS)  system,  which  played  a  major  role  in  the  CoBOP  field 
programs  (Davis  et  ah,  1999).  BRDFs  as  measured  by  K.  Voss  for  ooid  sand  and  grapestone  were 
used  for  simulations  of  bare  substrates,  and  various  analytical  BRDFs  developed  for  terrestrial  plant 
canopies  were  used  a  proxies  for  seagrass  BRDFs,  which  were  not  measured  in  the  CoBOP  field 
experiments. 

A  Backward  Monte  Carlo  model  was  developed  to  simulate  3-D  radiance  distributions  caused  by 
spatially  variable  bottom  reflectances  (such  as  patchy  sand  and  seagrass  bottoms)  and  by  sloping  or 
rippled  bottoms.  The  mathematical  techniques  are  described  in  Gordon  (1985),  Mobley  (1994, 

Section  6.2),  and  Mobley  and  Sundman  (2001c).  This  model,  called  BMC3D,  allows  the  user  to 
specify  an  arbitrary  spatial  pattern  of  bottom  BRDFs  (which  can  be  non-Lambertian)  as  well  as  sensor 
location,  viewing  direction,  and  field  of  view. 

Results  of  my  work  are  described  in  three  papers  that  have  been  submitted  to  the  CoBOP  special  issue 
of  Limnology  and  Oceanography  (Mobley,  Zhang,  and  Voss,  2001;  Mobley  and  Sundman,  2001c; 
Voss  et  ah,  2001).  In  addition  to  the  work  explicitly  described  here,  I  authored  or  co-authored  three 
other  papers  (Hoge  et  ah,  2000;  Leathers  et  ah,  2001;  Mobley,  Sundman,  and  Boss,  2001)  and  three 
book  chapters  (Bissett  et  ah,  2001;  Mobley  2001a,b)  that  contribute  to  the  overall  goals  of  CoBOP  or 
to  other  Navy  needs;  three  papers  were  presented  at  Ocean  Optics  XV. 

RESULTS 

Figure  1  (from  Mobley,  Zhang,  and  Voss,  2001)  shows  an  example  of  Hydrolight  simulations  of  the 
errors  that  occur  in  predictions  of  the  water-leaving  radiance  (the  quantity  of  interest  in  remote 

sensing)  if  a  non-Lambertian  bottom  is  replaced  by  a  Lambertian  bottom  having  the  same  irradiance 
reflectance  R.  The  non-Lambertian  bottoms  referenced  in  this  figure  include  a  measured  BRDF  of 
ooid  sand  at  the  CoBOP  field  site,  an  analytical  BRDF  with  parameters  chosen  to  simulate  sand 
(Hapke,  1993,  Eq.  8.89),  and  semianalytical  BRDFs  of  a  wheat  field  and  coniferous  forest  (Rahman,  et 
al,  1993),  which  were  used  as  proxies  for  a  seagrass  BRDF,  the  exact  form  of  which  is  unknown.  The 
water-column  optical  properties  were  taken  from  measured  values  at  the  CoBOP  Bahamas  site.  To 
accentuate  the  bottom  effects,  the  water  was  only  1  m  deep  and  the  sun  was  at  50  deg  from  the  zenith. 
These  errors  thus  represent  a  worse-case  situation  of  very  shallow  water  and  large  solar  angle.  The 
errors  become  smaller  as  the  bottom  becomes  deeper  or  the  sun  approaches  the  zenith.  0^,  and  (j)^  are 
the  polar  and  azimuthal  viewing  directions,  respectively.  The  viewing  directions  seen  in  a  PHILLS 
scan  are  shown  in  green.  The  sun  was  located  at  (0^,  4)^  =  (50°, 180°). 
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Fig.  1.  Contours  show  the  percent  errors  in  water-leaving  radiance  that  result  from 
replacing  a  non-Lambertian  BRDF  with  a  Lambertian  BRDF  having  the  same 
irradiance  reflectance.  Green  boxes  show  the  PHILLS  viewing  directions. 

[For  directions  relevant  to  remote  sensing,  the  errors  are  less  than  10%] 


Simulations  like  those  of  Fig.  1  yield  the  first  important  conclusion  of  my  work: 

For  solar  and  viewing  directions  relevant  to  most  ocean-color  remote  sensing,  a  non- 
Lambertian  bottom  can  be  modeled  as  a  Lambertian  bottom  having  the  same  irradiance 
reflectance  R  =  EJE^,  with  errors  of  less  than  10%  in  L^. 

Thus  the  different  angular  bottom  reflectance  patterns  do  not  make  much  difference  in  the  signals 
received  by  the  PHILLS  sensor.  The  crucial  measurement  that  optically  describes  a  bottom  type  is  the 
irradiance  reflectance. 

Figure  2  (from  Mobley  and  Sundman,  2001c)  shows  a  BMC3D  simulation  of  the  up  welling  irradiance 
E^,  at  the  surface,  as  would  result  from  moving  an  irradiance  sensor  across  a  boundary  between  bottom 
types  having  R  =  0.05  (e.g.,  seagrass)  and  R  =  0.5  (e.g.,  bright  sand).  E^  is  sometimes  used  in  remote 
sensing  studies  because  it  is  easily  measured  and  can  be  related  to  the  upwelling  radiance.  The  bottom 
was  at  a  depth  of  5  m,  and  the  water  lOPs  were  the  same  as  used  for  Fig.  1.  The  green  lines  in  the 
figure  show  the  Hydrolight-computed  values  for  infinite  grass  and  infinite  sand  bottoms  for 
comparison.  The  blue  line  shows  E^  as  computed  by  BMC3D.  In  this  simulation,  the  instrument 
begins  “seeing”  the  boundary  at  a  horizontal  distance  of  roughly  twice  the  water  depth.  The  red  line 
shows  E^  computed  using  the  simple  ID  model  of  Eqs.  (1)  and  (2)  below. 
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Fig.  2.  as  a  function  of  horizontal  position  near  a  boundary  between 
bottoms  of  different  reflectances. 


In  Fig.  2,  the  exact  (to  within  statistical  noise)  BMC3D  calculations  required  hours  of  computer  time; 
the  Hydrolight  calculations  took  less  than  5  seconds.  It  would  therefore  be  computationally  desirable 
to  have  an  approximate  model  that  could  use  the  Hydrolight  ID  solutions  for  the  two  different  bottoms 
in  combination  with  the  sensor  viewing  geometry  to  predict  (or  other  quantities)  near  such 
boundaries.  Such  a  ID  model  was  developed  from  simple  radiative  transfer  approximations.  The 
model,  when  applied  to  the  situation  seen  in  Figs.  2  and  3,  has  the  form 


Efx)  =  w{x-\)Ef\D-\)  +  w{x-2)EfW-2), 


(1) 


where 


/ 
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In  Eqs.  (1)  and  (2),  ^^(ID;/)  is  the  ID  (Hydrolight)  solution  for  bottom  type  i,  S(0)  is  the  sensor 
response  function  [S(0)  =  cos0  for  a  measurement  of  Ey],  c  is  the  beam  attenuation  of  the  water,  and 
the  integration  over  S(/)  means  “integrate  over  the  directions  where  the  sensor  at  location  x  sees 
bottom  type  As  seen  in  Fig.  2,  the  simple  ID  model  of  Eqs.  (1)  and  (2)  gives  rather  good 
predictions  of  the  3D  situation;  in  this  simulation,  the  difference  between  the  ID  model  and  the  exact 
3D  calculation  was  less  than  10%.  Additional  simulations  such  as  this  yield  the  second  important 
conclusion  of  this  work: 
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Using  a  ID  radiative  transfer  model  such  as  Hydrolight  to  predict  or  L„  can  give 
errors  of  many  tens  of  percent  near  the  boundary  of  highly  contrasting  bottom  types. 

However,  a  simple  analytical  model  that  weights  the  ID  solutions  according  to  the 
sensor  response  and  location  can  reduce  these  errors  to  roughly  10%. 

A  similar  statement  can  be  made  about  sloping  bottoms,  in  which  case  a  sloping  bottom  can  be 
modeled  as  a  level  bottom  using  Hydrolight  with  a  correction  for  the  change  in  the  incident  angle  of 
the  sun’s  unscattered  beam  onto  the  bottom. 

IMPACT/APPLICATION 

Inverse  models  are  being  developed  by  other  ONR  investigators  to  retrieve  bottom  depths  and  bottom 
classification  information  from  hyperspectral  ocean  color  sensors  such  as  PHILLS.  Those  retrieval 
algorithms  generally  assume  the  bottom  to  be  a  level,  Lambertian  reflector.  Radiative  transfer 
simulations  like  those  shown  here  allow  the  proposed  inversion  algorithms  to  be  evaluated  for  patchy, 
sloping,  non-Lambertian  bottom  boundaries. 

The  3-D  vs.  1-D  simulations  allow  a  determination  of  when  a  computationally  efficient  1-D  model  like 
Hydrolight  can  be  used  to  simulate  the  ocean  optical  environment,  and  when  a  much  more  expensive 
3-D  Monte  Carlo  model  must  be  run.  The  BMC3D  model  developed  here  also  allows  for  simulations 
of  particular  sensor  geometries  (such  as  sensor  angular  response),  which  cannot  be  exactly  simulated 
in  models  such  as  Hydrolight. 

TRANSITIONS 

1 .  The  BMC3D  code  was  transitioned  to  Dr.  Curtiss  Davis  of  NRL  Code  7212  for  use  in  computing  the 
effects  of  instrument  self-shading  on  the  buoyed  sensors,  which  are  used  to  acquire  sea-truth  data  for 
evaluation  of  hyperspectral  ocean-color  remote  sensing  systems.  Those  results  are  described  in 
Leathers  et  al.  (2001). 

2.  The  BMC3D  code  was  transitioned  to  Dr.  Vladimir  Haltrin  of  NRL  Code  7331  for  use  in  various 
simulations  of  underwater  3D  light  fields. 

3.  The  BMC3D  code  was  used  by  Mr.  Richard  Tompkins  of  the  Office  of  Naval  Intelligence  (ONI- 
241)  to  simulate  the  visibility  of  underwater  objects. 

4.  The  BMC3D  code  is  being  used  by  Todd  Bowers  at  the  Stennis  Space  Center  to  analyze  underwater 
visibility  data. 

RELATED  PROJECTS 

Measurements  of  BRDFs  for  various  bottom  types  were  made  by  K.  Voss  (Zhang  et  ah,  2001),  and 
measurements  of  water  and  bottom  optical  properties  and  associated  light  fields  were  made  by  various 
investigators  in  the  CoBOP  program.  Their  data  were  used  in  conjunction  with  my  models  to  evaluate 
the  results  of  the  optical  closure  experiment  conducted  during  the  2000  field  experiment  (Voss  et  ah, 
2001). 


5 


REFERENCES 


Bissett,  W.  P.,  O.  Schofield,  C.  D.  Mobley,  M.  F.  Crowley,  and  M.  A.  Moline,  2001.  Optical  remote 
sensing  techniques  in  biological  oceanography.  Invited  chapter  in  Methods  in  Microbiology  Vol.  30: 
Marine  Microbiology,  J.  H.  Paul,  Editor,  519-538. 

Davis,  C.  O.,  M.  Kappus,  J.  Bowles,  J.  Fisher,  J.  Antoniades,  and  M.  Carney,  1999.  Calibration, 
characterization  and  first  results  with  the  Ocean  PHILLS  hyperspectral  imager.  Proc.  SPIE  Conf.  on 
Imaging  Spectrometry  V,  SPIE  Vol.  3753:  160-168. 

Gordon,  H.  R.,  1985.  Ship  perturbation  of  irradiance  measurements  at  sea.  1:  Monte  Carlo 
simulations,  App/.  Opt.  24(23),  4172-4182. 

Hapke,  B.,  1993.  Theory  of  Reflectance  and  Emittance  Spectroscopy,  Cambridge  Univ.  Press,  455 
pages. 

Hoge,  F.  E.,  C.  D.  Mobley,  E.  K.  Sundman,  and  P.  E.  Eyon,  2000.  Radiative  transfer  equation 
inversion:  retrieval  of  oceanic  inherent  optical  properties.  J.  Geophys.  Res.,  submitted. 

Feathers,  R.  A.,  T.  V.  Downs,  and  C.  D.  Mobley,  2001.  Self-shading  correction  for  upwelling  sea- 
surface  radiance  measurements  made  with  buoyed  instruments.  Optics  Express,  8,  561-571. 

Mobley,  C.  D.,  1994.  Light  and  Water:  Radiative  Transfer  in  Natural  Waters.  Academic  Press,  592 
pages. 

Mobley,  C.  D.  2000a.  Radiative  Transfer  in  the  Ocean,  Invited  chapter  in  Encyclopedia  of  Ocean 
Sciences,  Academic  Press,  2001. 

Mobley,  C.  D.,  2001b.  Invited  contributions  to  Dictionary  of  Geophysics,  Astrophysics,  and 
Astronomy,  R.  Matzner,  editor.  CRC  Press. 

Mobley,  C.  D.  and  E.  K.  Sundman,  2001a.  Hydrolight  4.2  Users’  Guide.  Sequoia  Scientific,  Inc., 
Redmond,  WA,  87  pages. 

Mobley,  C.  D.  and  E.  K.  Sundman,  2001b.  Hydrolight  4.2  Technical  Documentation.  Sequoia 
Scientific,  Inc.,  Redmond,  WA,  76  pages. 

Mobley,  C.  D.  and  E.  K.  Sundman,  2001c.  Effects  of  optically  shallow  bottoms  on  upwelling 
radiances:  Effects  of  inhomogeneous  and  sloping  bottoms.  Limnol.  Oceanogr.,  submitted. 

Mobley,  C.  D.,  E.  K.  Sundman,  and  E.  Boss,  2001.  Phase  function  effects  on  oceanic  light  fields, 
Appl.  Optics,  submitted. 

Mobley,  C.  D.,  H.  Zhang,  and  K.  J.  Voss,  2001.  Effects  of  optically  shallow  bottoms  on  upwelling 
radiances:  Bidirectional  reflectance  distribution  function  effects.  Limnol.  Oceanogr.,  submitted. 


6 


Rahman,  H.,  B.  Pinty,  and  M.  M.  Verstraete,  1993.  Coupled  surface- atmosphere  reflectance  (CSAR) 
model  2.  Semiempirical  surface  model  usable  with  NOAA  AVHRR  data,  J.  Geophys.  Res.,  98(D11), 
20791-20801. 


Voss,  K.  J.,  C.  D.  Mobley,  L.  K.  Sundman,  J.  Ivey,  and  C.  Mazel,  2001.  The  spectral  upwelling 
radiance  distribution  in  optically  shallow  waters.  Limnol.  Oceanogr.,  submitted. 

Zhang,  H.,  K.  J.  Voss,  R.  P.  Reid,  and  E.  Louchard,  2001.  Bi-directional  reflectance  measurements  of 
sediments  in  the  vicinity  of  Lee  Stocking  Island,  Bahamas.  Limnol.  Oceanogr.,  submitted. 

PUBLICATIONS 

Max,  N.,  C.D.  Mobley,  B.  Keating,  and  E.-H.  Wu,  1997.  Plane-parallel  radiance  transport  for  global 
illumination  in  vegetation.  Rendering  ‘97,  J.  Dorsey  and  P.  Slusallek,  editors.  Springer  Verlag. 

Mobley,  C.D.  and  D.  Stramski,  1997.  Effects  of  microbial  particles  on  oceanic  optics:  Methodology 
for  radiative  transfer  modeling  and  example  simulations.  Limnol.  Oceanogr.,  42(3),  550-560. 

Stramski,  D.  and  C.D.  Mobley,  1997.  Effects  of  microbial  particles  on  oceanic  optics:  A  database  of 
single-particle  optical  properties.  Limnol.  Oceanogr.,  42(3),  538-549. 

Berwald,  J.,  D.  Stramski,  C.D.  Mobley,  and  D.A.  Kiefer,  1998.  The  effect  of  Raman  scattering  on  the 
average  cosine  and  diffuse  attenuation  coefficient  of  irradiance  in  the  ocean.  Limnol.  Oceanogr., 
43(4),  564-576. 

Lee,  Z.  P.,  K.  L.  Carder,  C.  D.  Mobley,  R.  G.  Steward,  and  J.  S.  Patch,  1998.  Hyperspectral  remote 
sensing  for  shallow  waters:  1.  A  semi-analytical  model.  Applied  Optics,  37(27),  6329-6338. 

Maffione,  R.A.,  J.M.  Voss,  and  C.D.  Mobley,  1998.  Theory  and  measurements  of  the  complete  beam 
spread  function  of  sea  ice.  Limnol.  Oceanogr.,  43(1),  29-33. 

Mobley,  C.D.,  G.E.  Cota,  T.C.  Grenfell,  R.A.  Maffione,  W.S.  Pegau,  D.K.  Perovich,  1998.  Modeling 
light  propagation  in  sea  ice.  IEEE  Trans.  Geosci.  Rem.  Sens.,  3(5(5),  1743-1749. 

Perovich,  D.K,,  J.  Longacre,  D.G.  Barber,  R.A.  Maffione,  G.E.  Cota,  C.D.  Mobley,  A.J.  Gow,  R.G. 
Onstott,  T.C.  Grenfell,  W.S.  Pegau,  M.  Landry,  and  C.S.  Roesler,  1998.  Eield  observations  of  the 
electromagnetic  properties  of  first-year  sea  ice,  IEEE  Trans.  Geosci.  Rem.  Sens.,  36(5),  1705-1715. 

Stephany,  S.,  E.  M.  Ramos,  H.  E.  de  Campos  Velho,  and  C.  D.  Mobley,  1998.  A  Methodology  for 
internal  light  source  estimation.  Computer  Model.  Simul.  Eng.,  3(3),  16 1-165. 

Lee,  Z.  P.,  K.  L.  Carder,  C.  D.  Mobley,  R.  G.  Steward,  and  J.  S.  Patch,  1999.  Hyperspectral  remote 
sensing  for  shallow  waters:  2.  Deriving  depths  and  optical  properties  by  optimization.  Applied  Optics, 
38(18),  3831-3843. 

Liu,  C-C,  J.  D  Woods,  and  C.  D.  Mobley,  1999.  Optical  model  for  use  in  oceanic  ecosystem  models, 
Appl.  Optics,  38(21),  4475-4485  . 


7 


Mobley,  C.  D.,  1999.  Estimation  of  the  remote-sensing  reflectance  from  above-surface  measurements. 
Appl.  Optics,  38(36),  7442-7455. 

Ohlmann,  J.  C.,  D.  A.  Siegel,  and  C.  D.  Mobley,  1999.  Ocean  radiant  heating:  1.  Optical  influences. 

J.  Phys.  Ocean.,  30,  1833-1848. 

Tyrrell,  T.,  P.M.  Holligan,  and  C.D.  Mobley,  1999.  Optical  impacts  of  oceanic  coccolithophore 
blooms.  J.  Geophys.  Res.,  104(C2),  3223-3241. 

Flatau,  P.  J.,  M.  Flatau,  J.  R.  V.  Zaneveld,  and  C.  D.  Mobley,  2000.  Remote  sensing  of  bubble  clouds 
in  seawater.  Quart.  J.  Royal  Meteor.  Soc.,  126(568),  2511-2524. 

Hoge,  F.  E.,  C.  D.  Mobley,  F.  K.  Sundman,  and  P.  E.  Fyon,  2000.  Radiative  transfer  equation 
inversion:  retrieval  of  oceanic  inherent  optical  properties.  J.  Geophys.  Res.,  submitted. 

Stephany,  S.,  F.  M.  Ramos,  H.  F.  de  Campos  Velho,  and  C.  D.  Mobley,  2000.  Identification  of 
inherent  optical  properties  and  bioluminescence  source  term  in  a  hydrologic  optics  problem.  J.  Quant. 
Spectros.  Rad.  Trans.,  67(2),  113-123. 

Stramska,  M.,  D.  Stramski,  B.  G.  Mitchell,  and  C.  D.  Mobley,  2000.  Estimation  of  the  absorption  and 
back-scattering  coefficients  from  in-water  radiometric  measurements.  Limnol.  Oceanog.,  45(3),  628- 
641. 

Bissett,  W.  P.,  O.  Schofield,  C.  D.  Mobley,  M.  F.  Crowley,  and  M.  A.  Moline,  2001.  Optical  remote 
sensing  techniques  in  biological  oceanography.  Invited  chapter  in  Methods  in  Microbiology  Vol.  30: 
Marine  Microbiology,  J.  H.  Paul,  Editor,  519-538. 

Feathers,  R.  A.,  T.  V.  Downs,  and  C.  D.  Mobley,  2001.  Self-shading  correction  for  upwelling  sea- 
surface  radiance  measurements  made  with  buoyed  instruments.  Optics  Express,  8,  561-571. 

Mobley,  C.  D.  2001.  Radiative  Transfer  in  the  Ocean,  Invited  chapter  in  Encyclopedia  of  Ocean 
Sciences,  Academic  Press,  2001. 

Mobley,  C.  D.,  2001.  Invited  contributions  to  Dictionary  of  Geophysics,  Astrophysics,  and  Astronomy, 
R.  Matzner,  editor.  CRC  Press. 

Mobley,  C.  D.  and  F.  K.  Sundman,  2001.  Effects  of  optically  shallow  bottoms  on  upwelling  radiances: 
Effects  of  inhomogeneous  and  sloping  bottoms.  Limnol.  Oceanogr.,  submitted. 

Mobley,  C.  D.,  H.  Zhang,  and  K.  J.  Voss,  2001.  Effects  of  optically  shallow  bottoms  on  upwelling 
radiances:  Bidirectional  reflectance  distribution  function  effects.  Limnol.  Oceanogr.,  submitted. 

Mobley,  C.  D.,  F.  K.  Sundman,  and  E.  Boss,  2001.  Phase  function  effects  on  oceanic  light  fields, 

Appl.  Optics,  submitted. 

Voss,  K.  J.,  C.  D.  Mobley,  F.  K.  Sundman,  J.  Ivey,  and  C.  Mazel,  2001.  The  spectral  upwelling 
radiance  distribution  in  optically  shallow  waters.  Limnol.  Oceanogr.,  submitted. 


8 


